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Abstract. 

A symplectic no-core shell model (Sp-NCSM) is constructed with the goal 
of extending the ab-initio NCSM to include strongly deformed higher- 
oscillator-shell configurations and to reach heavier nuclei that cannot be 
studied currently because the spaces encountered are too large to handle, 
even with the best of modern-day computers. This goal is achieved by 
integrating two powerful concepts: the ab-initio NCSM with that of the 
Sp(3, R) D SU(3) group-theoretical approach. The NCSM uses modern 
realistic nuclear interactions in model spaces that consists of many-body 
configurations up to a given number of hi} excitations together with mod- 
ern high-performance parallel computing techniques. The symplectic the- 
ory extends this picture by recognizing that when deformed configurations 
dominate, which they often do, the model space can be better selected so 
less relevant low-lying hQ configurations yield to more relevant high-lying 
h£l configurations, ones that respect a near symplectic symmetry found in 
the Hamiltonian. Results from an application of the Sp-NCSM to light 
nuclei are compared with those for the NCSM and with experiment. 



1 Introduction 

The concept of an ab initio no-core shell-model (NCSM) JTj, which yields a 
good description of the low-lying states in few-nucleon systems as well as in 
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more complex nuclei UlEl, nas taken center stage in the development of mi- 
croscopic tools for studying the structure of atomic nuclei. The architecture for 
the NCSM capitalizes on computational efficiencies that can be realized when 
many-particle Slater determinant basis states are mapped onto an integer bit 
string representation of that state on a computer. In addition, in the framework 
of the NCSM one can employ modern realistic interactions that reflect on the 
essence of the strong interaction. Recently developed realistic NN potentials 
include J-matrix inverse scattering potentials [3 1, high-precision NN potentials 
derived from meson exchange theory [4] and nuclear two- and many -body forces 
based on chiral effective field theory J5]. 

The symplectic no-core shell-model (Sp-NCSM) |!6 1 amplifies on this con- 
cept by recognizing that deformed configurations often dominate and these, 
while typically described by only few collective Sp(3,K) basis states, corre- 
spond to a special linear combination of a large number of NCSM basis states. 
Hence, the effective size of the model space can be significantly reduced and 
constrained to respect a near symplectic symmetry, that within a OhQ space 
reduces to SU(3), of the model Hamiltonian. In this way, the Sp-NCSM will 
allow one to account for even higher ft£l configurations required to realize ex- 
perimentally measured B(E2) values without an effective charge, and especially 
highly deformed spatial configurations required to reproduce a-cluster modes in 
heavier nuclei. 

As a 'proof-of-principle' study, results for no-core and symplectic no-core 
calculations up to 6M7 are compared for two nuclei, namely, the deformed 12 C 
and the closed-shell le O. The analysis of the results shows that the 0+ and the 
lowest 2 + and 4 + states in 12 C as well as the 0+ in le O, which are derived 
in the framework of the NCSM with the JISP16 realistic interaction U3J an d 
are well converged, reflect the presence of an underlying symplectic sp(3,R) 
algebraic structure^. This is achieved through the projection of realistic NCSM 
eigenstates onto Sp(3, M)-symmetric basis states of the symplectic shell model. 

The symplectic shell model Q |8] is a multiple oscillator shell general- 
ization of Elliott's SU(3) model and as well, a microscopic realization of the 
successful Bohr-Mottelson collective model. Symplectic algebraic approaches 
have achieved a very good reproduction of low-lying energies in 12 C using 
phenomenological interactions J9) or truncated symplectic basis with simplis- 
tic (semi-) microscopic interactions iflOl [TT1 . Here, we establish, for the first 
time, the dominance of the symplectic Sp(3, K) symmetry in nuclei as unveiled 
through ab initio calculations of the NCSM type with realistic interactions. This 
in turn opens up a new and exciting possibility for representing significant high- 
Ml collective modes by extending the NCSM basis space beyond its current 
limits through Sp(3,M) basis states, which yields a dramatically smaller basis 
space to achieve convergence of higher-lying collective modes. In this regard, 
it may be interesting to understand the importance of a larger model space be- 
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yond the 6hfl limit and its role in shaping other low-lying states in 12 C and le O 
such as the second + , which is likely to reflect a cluster-like behavior (e.g., see 
|[T2l ). This task, albeit challenging, is feasible for the no-core shell model with 
the symplectic Sp(3,R) extension. 

2 Symplectic Shell Model 

The symplectic shell model is based on the noncompact symplectic sp(3,R) 
algebra that with its subalgebraic structure unveils the underlying physics of a 
microscopic description of collective modes in nuclei OH). The latter follows 
from the fact that the mass quadrupole and monopole moment operators, the 
many-particle kinetic energy, the angular and vibrational momenta are all ele- 
ments of the sp(3,R) D su(3) D so(3) algebraic structure. Hence, collective 
states of a nucleus with well-developed quadrupole and monopole vibrations as 
well as collective rotations are described naturally in terms of irreducible repre- 
sentations (irreps) of Sp(3, R). Furthermore, the elements of the sp(3, R) alge- 
bra are constructed as bilinear products in the harmonic oscillator (HO) raising 
and lowering operators that in turn are expressed through particle coordinates 
and linear momenta. This means the basis states of a Sp(3, R) irrep can be ex- 
panded in a HO (m-scheme) basis, the same basis used in the NCSM, thereby 
facilitating symmetry identification. 

The symplectic basis states are labeled (in standard notation [7 8 1) according 
to the reduction chain 

Sp(3,R) D U(3) D SO(3) 

and are constructed by acting with polynomials V in the symplectic raising op- 
erator, A^ 20 \ on a set of basis states of the symplectic bandhead, |r CT }, which is 
a Sp(3, R) lowest-weight stated, 

\Y v Y n pY u K(LS)JMj)= [P T -{AW) x \Y a )] P ^ LS)JMj , (2) 

where T a = N a (X a labels Sp(3, R) irreps with (A ff fi a ) denoting a SU(3) 
lowest- weight state, r„ = n (A„ /i„), and = N w (A w p u ). The (A„ fj, n ) set 
gives the overall SU(3) symmetry of ^ coupled raising operators in V, (A w 
specifies the SU(3) symmetry of the symplectic state, and N u = N a + n is 
the total number of oscillator quanta related to the eigenvalue, N^hVt, of a HO 
Hamiltonian that is free of spurious modes. 

The symplectic raising operator A^ ', which is a SU(3) tensor with (\p) = 
(2 0) character, can be expressed as a bilinear product of the HO raising opera- 
tors, 



(211) 



(20) 

(3) 

lm 



'A Sp(3,R) lowest-weight state, |r CT ), is defined as J 4( 02 '|r CT ) = 0, where the symplectic 
lowering operator 2 ' is the adjoint of A^ 2 °' . 
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where the sums are over all A particles of the system. The first term in (0 de- 
scribes 2MI one-particle-one-hole (lp-lh) excitations (one particle raised by 
two shells) and the second term eliminates the spurious center-of-mass excita- 
tions in the construction (0. For the purpose of comparison to NCSM results, 
the basis states of the |r CT ) bandhead in (f2]i are constructed in a ?n-scheme basis, 

\T aK (L Q S Q )JoM Q ) = 

V^^\al)^^\al)\ X ^ |0), (4) 

where 0) is a vacuum state, V^ 71 ^ and Vg X " ^ denote polynomials of pro- 
ton (ap and neutron (aj,) creation operators coupled to good SU(3)xSU(2) 
symmetry. 

3 Results and Discussions 

The lowest-lying states of 12 C and le O were calculated using the NCSM as im- 
plemented through the Many Fermion Dynamics (MFD) code iTPTl . For 12 C 
we used an effective interaction derived from the realistic J1SP16 AW potential 
(0 for different W. oscillator strengths, while for le O the bare JISP16 interac- 
tion was employed. We are particularly interested in the J = 0+ and the lowest 
J = 2 + (= 2^) and J = 4+ (=4+) states of the ground-state (gs) rotational band 
in 12 C and the J = 0+ s state in le O that appear to be well converged in the 
Nmax = 6 NCSM basis space. 

Here we report on an analysis that is restricted to Op-Oh configurations. It is 
important to note that 2ft£l 2p-2h (2 particles raised by one shell each) and higher 
rank np-nh excitations and allowed multiples thereof can be included by build- 
ing them into an expanded set of lowest- weight Sp(3, R) starting state configura- 
tions. The same "build-up" logic, (ffj), holds because by construction these addi- 
tional starting state configurations are also required to be lowest- weight Sp(3, R) 
states. Note that if one were to include all possible lowest-weight np-nh start- 
ing state configurations (n < N max ), and allowed multiples thereof, one would 
span the entire NCSM space. The addition of 2WI 2p-2h, AW 4p-4h, and 
higher configurations, which build upon more complex starting states, will be 
the subject of a follow-on investigation. 

3.1 Ground-state rotational band in the 12 C nucleus 

For 12 C there are 13 unique Op-Oh Sp(3,R) irreps which form the symplectic 
bandhead basis states, \T a ) with N a = 24.5. For each Op-Oh Sp(3, R) irrep we 
generated basis states according to (ffji up to N max = 6 (6fi.il model space). 
The typical dimension of a symplectic irrep basis in the N max = 6 space is on 
the order of 10 2 as compared to 10 7 for the full NCSM m-scheme basis space. 

As N max is increased the dimension of the J = 0,2, and 4 symplectic 
space built on the Op-Oh Sp(3, R) irreps grows very slowly compared to the 
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NCSM space dimension (Fig. Q]). This means that a space spanned by a set of 
symplectic basis states may be computationally manageable even when high- 
fill configurations are included. 

Analysis of overlaps of the symplectic states with the NCSM eigenstates for 
the 0+ and the lowest 2 + and 4 + states reveals nonnegligible overlaps for only 
3 of the 13 Op-Oh Sp(3,R) (N a = 24.5) irreps, specifically, the leading (most 
deformed) representation (A CT /v) = (0 4) carrying spin S = together with 
two 5 = 1 representations with identical labels (12) but different bandhead 
constructions for protons and neutrons namely, {(A„- /Lt 7r )5 7r , {\ v fJ, v )S v } is 
{(0 2)0, (1 0)1} and {(1 0)1, (0 2)0}. The dominance of only three irreps 
additionally reduces the dimensionality of the symplectic model space (Fig. [T] 
red diamonds). 

The overlaps of the most dominant symplectic states with the NCSM eigen- 
states for the 0+,, 2^~ and 4+ states in the 0, 2, 4 and 6hil subspaces are given 
in Table [T] The results show that approximately 80% of the NCSM eigenstates 
fall within a subspace spanned by the 3 leading Op-Oh Sp(3, R) irreps, with the 
most deformed irrep, (0 4), carrying about 65% of the 80%. In order to speed 
up the calculations, we retained only the largest amplitudes of the NCSM states, 
those sufficient to account for at least 98% of the norm which is quoted also in 
the table. 

In addition, the 0+ analysis of the 5 = (5 = 1) part of the NCSM wave- 
function reveals that within each Ml subspace only about 1 — 1.5% of the NCSM 
0+. are not accounted for by the 5 = (5 = 1) Sp(3, R) irrep(s) under con- 
sideration. In the N max = 6 model space the 5 = symplectic irrep and the 
two 5=1 irreps account for 91% and 80%, respectively, of the corresponding 
5 = and 5=1 parts of the NCSM realistic eigenstate for the J = 0+ in 12 C. 
In summary, the 5 = plus 5=1 part of the NCSM wavefunction is very well 




Figure 1. Dimension of the NCSM (blue squares) and J — 0,2, and 4 Sp(3, R) (red 
diamonds for the 3 most significant Op-Oh irrep case and green circles for when all 13 Op- 
Oh irreps are included) model spaces as a function of maximum allowed ftfi excitations, 

max ■ 
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explained by only the three Sp(3, R) collective configurations. 

How the results presented in Table Q] change as a function of the oscillator 
strength Ml is shown in Fig. |2]for the case of the 0+ s state. Clearly, the projec- 
tion of the NCSM wavefunctions onto the symplectic space slightly changes as 
one varies the oscillator strength Ml. The 3 Sp(3, R) irreps, (0 4)5 = and the 
two (1 2)5 = 1, remain dominant, only their contributions change. The overall 
overlaps increase towards smaller Ml HO frequencies and, for example, for 0+ s 
it is 85% in the N max = 6 and Ml = HMeV case. Clearly, the largest con- 
tribution comes from the leading, most deformed, (0 4)5 = Sp(3,R) irrep, 
growing to 91% of the total Sp(3, R)-symmetric part for Ml =1 1 MeV. As ex- 
pected, Fig. |2] also confirms that with increasing Ml the higher Ml excitations 
contribute less while the lower QMl configurations grow in importance. 

In short, the low-lying states in 12 C are well described in terms of only three 
Sp(3,R) irreps with total dimensionality of 514, which is only 0.001% of the 
NCSM space, with a clear dominance of the most deformed (0 4)5 = col- 
lective configuration. It is important to note that our results suggest that over- 
laps can be further improved by the inclusion of the most important 2MI 2p-2h 
Sp(3, R) irreps. In this way it may be possible to achieve overlaps of more then 
90% while keeping the size of the basis space small, possibly much less than 
1 % of the NCSM result. This is the subject of ongoing investigations and will 
be addressed in a subsequent study. 



Table 1. Probability distribution of NCSM eigenstates for C across the leading 3 Op-Oh 
Sp(3, R) irreps, hn=\5 MeV. 
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Figure 2. 12 C ground + state probability distribution over 0M7 (blue, lowest) to 
6h£l (green, highest) subspaces for the 3 Op-Oh Sp(3, R) irrep case (left) and NCSM 
(right) together with the (0 4) irrep contribution (black diamonds) as a function of the 
hSl oscillator strength in MeV for N max = 6. 



The 0+,, 2^ and 4^ states, constructed in terms of the three Sp(3, R) irreps 
with probability amplitudes defined by the overlaps with the NCSM wavefunc- 
tions, were also used to determine B(E2) transition rates. The B(E2 : 2^ — > 
0+ s ) value, for example, turns out to be «110% of the corresponding NCSM 
number for the Ml = 15MeV and N max = 4 case. While this ratio decreases 
slightly for smaller Ml oscillator strengths, it is significant that this estimate 
for the dominant Sp(3, R) configurations exceeds the corresponding full NCSM 
results and therefore lies closer to the experimental B(E2 : 2^ — > 0+ ) value. 

3.2 Ground state in the 16 nucleus 

The Sp-NCSM is also applied to the ground state of a closed-shell nucleus 
like le O. There is only one Op-Oh Sp(3,R) irrep with spin 5 = and T a 
specified by N a = 34.5 and (A CT n„) = (0 0). As in the 12 C case, for the 
Op-Oh Sp(3, R) irrep we generated basis states according to (O up to N max = 6 
(6MI model space), which yields a symplectic model space that is only a 
fraction (w 0.1%) of the size of the NCSM space. Consistent with the outcome 
for 12 C, the projection of the NCSM eigenstates onto the symplectic basis 
reveals a large Sp(3, R)-symmetric content in the ground-state wavefunction 
(Fig. |3). Furthermore, the overall overlap increases by w 10% when the most 
significant 2MI 2p-2h are included. 

While the focus here has been on demonstrating the existence of Sp(3, R) 
symmetry in NCSM results for 12 C and le O, and therefore a possible path for- 
ward for extending the NCSM to a Sp-NCSM scheme, the results can also be 
interpreted as a further strong confirmation of Elliott's SU(3) model since the 
projection of the NCSM states onto the Qhil space [Fig. |2]and Fig[3] blue (right) 
bars] is a projection of the NCSM results onto the SU(3) shell model. For le O 
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Sp(3,R) NCSM 



Figure 3. O ground state probability distribution over OhQ (blue, lowest) to 
QhQ. (green, highest) subspaces for the leading Op-Oh (0 0) Sp(3, R) irrep case (left) 
and NCSM (right) for N max = 6 and bare JISP16 interaction. 

the Ohfl SU(3) symmetry is w 60% of the NCSM 0+ [Fig. [3] blue (left) bars]. 
For 12 C the 0M2 SU(3) symmetry ranges from just over 40% of the NCSM 0+ 
for Ml = 1 1 MeV to nearly 65% for Ml =18 MeV [Fig. El blue (left) bars] with 
80%-90% of this symmetry governed by the leading (04) irrep. These num- 
bers are consistent with what has been shown to be a dominance of the leading 
SU(3) symmetry for SU(3)-based shell-model studies with realistic interactions 
in QMl model spaces. It seems the simplest of Elliott's collective states can be 
regarded as a good first-order approximation in the presence of realistic interac- 
tions, whether the latter is restricted to a QMl model space or the richer multi-M7 
NCSM model spaces. 

4 Conclusions 

Wavefunctions, which are obtained in a NCSM analysis with the JISP16 realistic 
interaction, project at approximately the 80% level onto the leading (three) Op- 
Oh irreps of the corresponding Sp-NCSM for the lowest 0+,, 2f and 4+ states 
in 12 C and at more than 70% level for the ground state in the closed-shell le O 
nucleus. (While not part of the current analysis, preliminary results indicate that 
when the space is expanded to include the most important 2MI 2p-2h irreps the 
percentage grows by approximately 10%.) The results confirm for the first time 
the validity of the Sp(3, M.) approach when realistic interactions are invoked and 
hence demonstrate the importance of the Sp(3, R) symmetry in light nuclei as 
well as reaffirm the value of the simpler SU(3) model upon which it is based. 

The results further suggest that a Sp-NCSM extension of the NCSM may be 
a practical scheme for achieving convergence to measured B(E2) values without 
the need for introducing an effective charge and even for modeling cluster-like 
phenomena as these modes can be accommodated within the general framework 
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of the Sp(3, R) model if extended to large model spaces (high N max ), but with 
a size that is typically only a fraction of the NCSM size. This suggests that a Sp- 
NCSM code could allow one to extend no-core calculations to higher Ml con- 
figurations and heavier nuclei that are currently unreachable because the model 
space is typically too large to handle, even on the best of modern day computers. 
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